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Abstract. Using a new, increased dataset of 7 QSOs from VLT/UVES observations combined with one QSO 
from the literature, the minimum Doppler parameters as a function of neutral hydrogen column density Nhi, 
b c {Nni), of the Lya forest has been derived at three redshifts <z> — 2.1, 3.3 and 3.8. In particular, five QSOs at 
<2>= 2.1 enable us to study the cosmic variance of ^(A^hi) at lower z for the first time. When incompleteness 
of the number of the observed lines towards lower Ahi is accounted for, the derived slopes of b c (Nni), (T — 1), are 
consistent with no- z evolution with an indication of lower value at < z > = 3.3, while & c (Ahi) at a fixed column 
density Am = 10 136 cm" 2 , 6 C (13.6), increases as z decreases. Assuming a QSO-dominated UV background, the 
slope of the equation of state (7 — 1) shows no z-evolution within large uncertainties and the temperature at 
the mean density, To, decreases as z decreases at three redshift ranges. There is a large fluctuation of (T — 1) 
and & c (13.6) even at the similar redshifts, in particular at < z > — 3.3 and 3.8. The lower (r — 1) and higher 
6 C (13.6) values at z ~ 3.1 and 3.6 compared to ones at z ~ 3.4 and 3.9 are caused by a lack of lower- ATm and 
lower-6 lines at lower- z parts of each QSO at z > 3, probably due to the Hell reionization. This result suggests 
that an impact from the Hell reionization on the forest might be mainly on the lower- Ah 1 forest. From this new 
dataset, we find some forest clouds with a high ratio of Si IV column density to C IV column density, Nsiw/Nc iv, at 
z < 2.5, although the bulk of the forest clouds shows lower Nsiw/Ncw- This high Asiiv/Afciv at z < 2.5 suggests 
that some forest clouds are exposed to a soft UV background. This lack of strong discontinuity of Nsiw / New at 
A^hi = 10 14_17 cm -2 at z ~ 3 suggests that Asiiv/Afciv might not be a good observational tool to probe the Hell 
reionization and/or that the UV background might be strongly affected by local, high-z galaxies at z < 3. 
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1. Introduction 

The Lya forest imprinted in the spectra of high-z 
QSOs arises from the fluctuating low-density intergalactic 
medium (IGM), highly photoionized by the metagalactic 
UV background. Since the universe expands adiabatically 
and the Lya forest is in photoionization equilibrium with 
the UV background, the temperature of the Lya forest 
as a function of z provides a unique and powerful tool to 
probe the physical state of the IGM and the r eionization 
history of the universe (Hui & Gnedin 1997 ; Schaye et 
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* The data used in this study are based on public data 
released from the VLT/UVES Commissioning and Science 
Verification and from the OPC program 65.O-0296A (P.I. 
S. D'Odorico) at the VLT/Kueyen telescope, ESO, Paranal, 
Chile. 
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For a low-density (the baryon overdensity S < ~ 10), 
photoionized gas, the temperature of the gas is shown to 
be tightly correlated with the overdensity of the gas. This 
relation, i.e. the equation of state, is defined by T = T (l + 
(5) 7_1 , where T is the gas temperature in K, T is the gas 
temperature in K at the mean gas density and (7 — 1) is 
a constant at a given redshift z. Both To and (7 — 1) are 
a function of z, depending on the thermal history of the 
IGM (Hui & Gnedin |1997|) . 

This equation of state, however, is not directly observ- 
able. Instead of T and (1 + 6), observations only provide 
the neutral hydrogen column density A/ Hl (in cm" 2 ) and 
the Doppler parameter b (in km s" 1 ) of the forest ab- 
sorption lines. In practice, a lower cutoff envelope in the 
Nn r~b distribution is used to probe the upper limit on the 
temperature of the IGM since the forest lines could be 
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observations and data reduction, and the fitted line lists 



may be found in Kim et al. ( |200Ia| , |2001b| ). 

The line parameters of Q0000-263 were taken from 



broadened by processes other than the thermal broaden- 
ing. Translating a -/Vhi-6 envelope into a (I + S)-T rela- 
tion depends on many physical assumptions, such as the 
ionizing UV backg round J„ (Miralda-Escude et al. 1996 ; Lu et al. (p.996|) to include the highest redshift Lya forest 
Schaye et al. |1999| ). 

This minimum Doppler cutoff 6 c (^hi) can be de- 
scribed by 



log(& c ) = bg(6o) + (T - 1) log(JV Hl ), 



(1) 



where log(&o) is the intercept of the cutoff in the log(ATm)- 
log b diagra m and (T — 1) is the slope of the cutoff (Schaye 
et al. (1999b . 

From observations alone, both n o z-ev olution of iVm- 
independent b c (Kirkman & Tytler 1997 ; Savaglio et al. 
1999] ) and increasing b c with decreasing z (Kim et al. 
199JJ) have been claimed. Results from simulations com- 
bined with observations have also claimed both no-z evo- 
lution of T and (7 - 1) (McD onald et al. |2000|) and a 



z-evolution (Ricotti et al. 2000 ; Schaye et al. [20001 Kim, 
Cristiani & D'Odorico 2001a ). Deriving 6 c (A^ji) from ob- 



servations depends on many factors such as the method of 
line deblending, the number of available absorption lines, 
the mctal-linc contamination, and the method of fitting 
the lower iVm-& envelope (Hu e t al. |1995 ; Kirkman & 
Tytle r |1997| Bryan fc M achacek ^OOo[ McD onald et al. 
200Cj; Ricotti et al. 



200C 



200C; Kim et al. 



Shaye et al. 

2001a). The different approaches and the limited numbers 
of lines have led, in part, to the contradicting results on 
the evolution of b c (Nni) in the literature. 

Here, using a new, increased dataset from 7 QSOs ob- 
served with the VLT/UVES combined with the published 
data on one QSO obtained with Keck/HIRES, we present 
the evolution of the Doppler cutoff 6 c (-/Vhi) at three red- 
shifts < z > = 2.1, 3.3 and 3.8. In particular, five QSOs 
at < z > = 2.1 enable us to study the cosmic variance 
of 6 c (A^hi) and to improve a determination of b c (N-m) at 
lower z for the first time. In Sect. 2, we briefly describe the 
data used in this study. The analyses of the observations 
are presented in Sect. 3. The discussion is in Sect. 4 and 
the conclusions are summarized in Sect. 5. In this study, 
all the quoted uncertainties are la errors. 

2. Data used in the study 

Table |l| lists the eight QSOs analyzed in this paper. 
Seven QSOs were observed with the UVES spectro- 
graph at the VLT Kueyen telescope (built by ESO, P.I. 
S. D'Odorico). The UVES data were reduced with the 
MIDAS ECHELLE/UVES package. The final reduced 
vacuum heliocentric spectra have S/N of 30-50 per pixel 
in the regions of interest and a resolution of R ~ 45, 000. 
The spectra were normalized locally using a 5th order 
polynomial fit. The normalized spectra were then fit- 
tcd with Voigt profiles using VPFIT (C arswell et al.: 
http: / /www. ast.cam.ac.uk/~rfc/vpfit. html ) with the re- 
duced x 2 threshold of 1.3 to obtain the three line param- 
eters, z, b and JVhi- The metal lin es were identified and 
removed as described in Kim et al. ( 2001a ). Details of the 



available in the literature, with similar resolution and S/N 
to the UVES data. Their analysis of Q0000-263 was also 
undertaken with VPFIT. 

In order to avoid confusion with the Ly/3 forest and the 
proximity effect, we consider only the wavelength range 
from the Ly/3 emission to 3 000 km s _1 shortward of the 
Lya emission. However, the redshift intervals actually used 
are further limited by other factors such as the incomplete 
coverage of the forest region, a damped Lya system and 
our attempt to overlap the wavelengths of each QSO as 
much as possible to study the cosmic variance of 6 c (-/Vh~i). 
Table [l] lists the wavelength ranges used for each QSO. 

We restrict our analysis to iVm = io 12 - 5-14 - 5 cm~ 2 . 
The lower limit corresponds to the detection threshold in 
the regions of poorest S/N and the upper limit is where 
the A^hi estimate from fitting Lya alone becomes unreli- 
able because of line saturation. Because lines in blends can 
also have large uncertainties, we have further restricted 
the analysis to include only those lines with profile fit- 
ting errors less than 25% in A^hi and b to be tter define 
the low er cutoff envelopes (Schaye et al. 2000 
|2001aD . 



Kim et al. 



In this study, Sample A defines all the lines available 
from all QSOs which have treated as a single dataset at 
each z. We also define Sample B in order to study a fluctu- 
ation of the Doppler cutoff at similar redshifts. The spec- 
tral coverage for each QSO from the same z bin is not 
uniform. For those QSOs with more than ~ 600A cover- 
age (HE1122-1648, HE2217-2818, Q0055-269 and Q0000- 
263), the line lists have been divided into two subsets: a 
group at higher redshifts and a group at lower redshifts. 
The rest of the QSOs do not have enough coverage to 
make this splitting possible and provide only one group 
each. We label the ensemble of these groups Sample B. 
Each group of Sample B spans 300A-350A and is defined 
to have roughly a similar redshift coverage. 



3. The physical state of the IGM 

3.1. The z-evolution of the Doppler cutoff 

To derive 6 c (Ahi), we adopted two methods: the iterative 
power-law fit (Schaye et al. 1999, 2000|) a nd the smoothed 
b power-law fit (Bryan & Machacek 2000] ) . From these two 
power-law fit methods, we derived log(&o,fc) and (T— l)k of 
Eq. |l| for k — i or s, where subscripts i and s indicate the 
iterative power-law fit and the smoothed b power-law fit, 
respectively. For the iterative power-law fit, 100 bootstrap 
realizations were averaged with an exclusion threshold of 
0.5 absolute mean deviation (AMD)[]. For the smoothed 



b power-law fit, a smoothing constant of 3 km s was 



1 Schaye et al. (2000) adopted a threshold of 1 AMD, which 
results in a power law fit with higher \ 2 ■ 
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Table 1. Analyzed QSOs 



QSO 


2 a 
"cm 


mag a 


AA (A) 




ya 


# of lines b 


Comments 


Q1101-264 


2.145 


16.0 


3500-3778 


1.88- 


-2.11 


69 


UVES SV, a damped system atz= 1.8386 


J2233-606 


2.238 


17.5 


3500-3890 


1.88- 


-2.20 


88 


UVES Commissioning I 


HE1122-1648 


2.400 


17.7 


3500-4091 


1.88- 


-2.37 


179 


UVES SV, split into 2 C 


HE2217-2818 


2.413 


16.0 


3510-4100 


1.89- 


-2.37 


159 


UVES Commissioning I, split into 2 d 


HE1347-2457 


2.534 


16.8 


3760-4100 


2.09- 


-2.37 


91 


UVES SV, incomplete observations 


Q0302-003 


3.281 


18.4 


4808-5150 


2.96- 


-3.24 


107 


UVES Commissioning I, incomplete observations 


Q0055-269 


3.655 


17.9 


4850-5598 


2.99- 


-3.60 


264 


UVES, Sept 20-22, 2000, split into 2 C 


Q0000-263 


4.127 


17.9 


5450-6100 


3.48- 


-4.02 


209 


Lu et al. (199£ ), split into 2 f 



a Taken from the SIMBAD database. The magnitude of HE1347-2457 is from NED. 
b For Nm = 10 12 5 " 14 - 5 cm~ 2 . Only for lines with the errors less than 25% in both AT H , and b. 
c For Sample B, the spectrum is split into 3500-3800 A (87 lines) and 3800-4091 A (92 lines). 
d For Sample B, the spectrum is split into 3510-3800 A (76 lines) and 3800-4100 A (83 lines). 
c For Sample B, the spectrum is split into 4850-5220 A (110 lines) and 5220-5598 A (154 lines) 
f For Sample B, the spectrum is split into 5450-5820 A (106 lines) and 5820-6100 A (103 lines) 



Table 2. The power-law fits to b c (Nnj 



Sample A 



<z> 


log 




# of lines 


log(feo.i) 


(r-i)i 


6o,i,(13.6) 


log(6 , s ) 


(r - l)s 


b c ,s(13.6) 


2.1 


13.0- 


-14.5 


349 


-0.745 ± 0.089 


0.150 ±0.006 


19.8 ±0.8 


-0.495 ± 0.074 


0.131 ±0.002 


19.1 ± 1.1 


3.3 


13.0- 


-14.5 


275 


-0.413 ±0.116 


0.122 ±0.008 


17.2 ± 1.0 


0.279 ± 0.473 


0.072 ±0.127 


18.4 ± 1.5 


3.8 


13.3 


-14.5 


152 


-0.948 ±0.125 


0.159 ±0.008 


16.8 ± 1.6 


-2.699 ±0.159 


0.285 ±0.042 


14.9 ± 1.1 


2.0 a 


13.0- 


-14.5 


176 


-0.214 ±0.180 


0.111 ±0.013 


19.1 ± 1.1 


-1.078 ±0.173 


0.172 ±0.047 


18.0 ± 1.1 


2.2 a 


13.0- 


-14.5 


173 


-0.954 ±0.132 


0.167 ±0.010 


19.7 ± 1.0 


0.400 ± 0.110 


0.066 ± 0.030 


19.9 ± 1.1 


3.1 b 


13.0- 


-14.5 


157 


0.645 ± 0.379 


0.048 ± 0.027 


19.9 ± 1.5 


-0.366 ± 0.103 


0.122 ±0.028 


19.6 ± 1.1 


3.4 b 


13.0- 


-14.5 


118 


-0.330 ± 0.053 


0.117 ±0.009 


17.5 ±0.8 


-0.745 ±0.131 


0.142 ±0.035 


15.5 ± 1.1 


3.6 C 


13.3 


-14.5 


74 


0.572 ± 0.331 


0.054 ± 0.023 


20.4 ± 0.8 


1.092 ± 0.007 


0.014 ±0.002 


19.0 ± 1.0 


3.9 C 


13.3 


-14.5 


78 


-0.697 ±0.133 


0.143 ±0.010 


17.5 ±0.7 


-0.499 ± 0.235 


0.122 ±0.062 


14.5 ± 1.1 



Sample B (averaged for the individual members) 



2.1 

3.3 
3.8 



13.0-14.5 
13.0-14.5 
13.3-14.5 



-0.044 ± 0.506 
0.617 ±0.966 
-0.063 ± 0.898 



0.102 ±0.037 
0.051 ± 0.067 
0.098 ± 0.062 



21.4 ±0.8 
20.3 ±2.7 
19.0 ± 2.1 



0.053 ± 0.731 
0.283 ± 0.970 
0.182 ± 1.350 



0.090 ± 0.053 
0.072 ± 0.068 
0.078 ± 0.094 



19.2 ± 1.1 

18.3 ±2.7 

17.4 ±3.1 



Results from Schaye et al. (200C) (their sample corresponds to our Sample B) 



3.1 
3.1 
3.8 



12.5-14.5 
12.5-14.5 
12.5-14.8 



< 0.150 
~ 

< 0.150 



23-24 
20-22 
18-19 



In order to have a similar number of lines at z ~ 2.0 and z ~ 2.2, the < z > = 2 sample consists of the line parameters 
from Q0011-264, J2233-606, HE1122-1648 at 3500-3760 A and HE2217-2818 at 3510-3760 A. The the <z> = 2.2 sample 
consists of the line parameters from Q1347-2457, HE1122-1648 at 3760-4100 A and HE2217-2818 at 3760-4100 A. 
The < z > = 3.1 sample and the < z > = 3.4 sample are taken from Q0302-003 and Q0055-269 at 4850-5220 A, from 
Q0055-269 at 5220-5598 A, respectively. 

The < z > = 3.6 sample and the <z> — 3.9 sample are taken from 5450-5820 A and from 5820-6100 A, respectively. 



used for each subsample having 30 lines and the robust 
power-law fit was applied. 

Schaye et al. ( [2000 ) used the iterative power-law fit 
for iV H i = 10 12 - 5 ^^cm- 2 at 2 < z < 3.7 and for 
JVhi = 10 12 - 5 - 14 - 8 cm- 2 at 3.7 < z < 4.4. Lines with 
^Yhi < 10 13 cm~ 2 , however, suffers from incompleteness 
due to line blending. For example, at z ~ 2.1, lines with 
^Yhi < 10 13 cm -2 are affected by incompleteness and this 



threshold increases with z (cf. Hu et al. 1995| ; Kim et al 



1997, 2001a, 2001t). Incompleteness causes a bias in the 
measurements of (T — 1). Due to the lack of lines with 
Vhi below this threshold and b < 15 km s" 1 (cf. Kim 
et al. 2001a), an observationally measured (r — 1) value 
becomes natter than its true value. We have defined dif- 
ferent lower TVh i fitting thresholds at the various redshifts 
in order to avoid this bias (see Fig. and to obtain a sta- 
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Fig. 1. The Ahi~& diagrams at < z > = 2.1, 3.3 and 3.8. Errors are not displayed. Triangles indicate possible metal 
lines or lines in the Hi complex mixed with metal lines. Dotted lines indicate the Ahi ranges considered in the 
study. Dot-dashed lines represent the lower Ahi fitting threshold actually used in the fit, above which incompleteness 
is negligible. Solid lines and dashed lines represent the iterative power-law fits and the smoothed b power-law fits, 
respectively. Shaded area represents a N-Rj-b distribution enclosed by two power-law fits at <z> = 2.1. 



ble estimate of (T — 1) and b c at the fixed column density 
N ni = 10 13 - 6 cm- 2 , 6 C (13.6). 

Fig. [l] shows the Ahi-6 diagram at < z > = 2.1, 3.3 
and 3.8 for Sample A. The fitted parameters are listed 
in Table as well as the Ahi ranges used in the fit. At 
higher z, several lines show a b value smaller than 6 c (Ahi) 
at < z >— 2.1 (shaded area). In fact, 6 C (13.6) from the 
both power-law fits increases as z decreases. The slopes 
(T — l)i at < z > = 2.1 and 3.8 are similar within errors, 
while the slope (r — l) s at < z > = 2.1 is flatter than 
the one at < z > — 3.8 more than 3a. Both slopes shows 
the lowest value at < z > = 3.3 (see Table ||). Note that, 
however, the errors from the bootstrap method are likely 
to be underestimated (cf. Schaye et al. 2000). The close 
examination of Fig. |l| suggests that (r — l) s at < z > = 3.8 
might well be overestimated and that the real (r— 1) might 
be inbetween (r — 1); and (r — l) s . 

We visually compared our 6 c j(13.6) with the results of 
Schaye et al. ( 2000| ; their figure 1), although their sample 
corresponds to our Sample B. As shown in Table g, their 
6 C (13.6) values are larger than our 6 C (13.6) values at simi- 
lar z with the differences larger at smaller z. Their (r — 1); 
values appear to be flatter than our values at all z, likely 
caused by their lack of accounting for incompleteness of 
observed lines. 

3.2. Fluctuations of the Doppler cutoff 

Fig. | shows (r - l)j and 6 C)i (13.6) as a function of the 
lower iVm threshold iVm th f° r Sample A. The lack of 
lines due to incompleteness results in flatter (r — l)i at 
lower iVHi.th- When the lower JVm threshold increases so 
that incompleteness does not affect the Ah i-b distribution 
severely, (r — l)j becomes stabilized (the lower JVm,th ~ 
10 13 cm -2 at <z> — 2.1 and 3.3, and the lower iVni,th ~ 
10 13 - 3 cm -2 at < z > = 3.8). If the available lines in the 
fit, however, become too small for higher A^i.thi (r — l)i 
becomes rather ill-defined. On the other hand, the lower 
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Fig. 2. Slopes (r — l)j and 6 c ,i(13.6) as a function of a 
lower Ah i threshold Am^h for Sample A. Open circles, 
open squares and open diamonds indicate < z > = 2.1, 3.3 
and 3.8, respectively. 



Am.th does not affect the & Ci i(13.6) values since 6 Cj ; at 
Ah i ~ 10 13 6 cm -2 behaves more like a pivotal point. The 
slopes (r — l) s and 6 CjS (13.6) from the smoothed b power- 
law fit also show similar behaviors. 

Fig. H shows (r — l)j and 6 Cj i(13.6) as a function of 
z for the individual members of Sample B (open sym- 
bols), the averaged values from Sample B (filled triangles) 
and Sample A (filled symbols). The members of Sample B 
show the large fluctuation of (r — 1), and & c .i(13.6) even 
at the similar z (see Table ^). With more lines available 
to sharpen the lower cutoff envelope, (r — l)i becomes 
steeper than the ones derived from the individual mem- 
bers using a smaller number of lines, possibly approaching 
asymptotic values (Sample A: filled circles, filled squares 
and filled diamonds). For both Sample A and Sample B 
averaged (filled triangles), (r — l)i shows the lowest value 
at < z > = 3.3, although uncertainties for Sample B av- 
eraged are rather large at <z> — 3.3 and 3.8. While the 
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Fig. 3. Slopes (T — l)j and 6 c ,i(13.6) as a function of z. 
Open symbols (circles, squares and diamonds at < z > = 
2.1, 3.3 and 3.8, respectively), filled symbols (same as in 
open symbols) and filled triangles represent the individual 
members of Sample B, Sample A and Sample B averaged 
at each z, respectively. 



lcr is 0.037 at < z > = 2.1, it becomes twice as large as 
that at < z > — 3.3 and 3.8. This larger fluctuation at 
z > 3 is caused in part by a smaller data size and in part 
by the cosmic variance. The fluctuation at < z > = 2.1 is 
mostly caused by different sightlines. A single, long sight- 
line does not show any significant difference between the 
lower- z part and the higher-z part, such as HE1122-1648 
and HE2217-2818. On the other hand, at < z > = 3.3 
and 3.8, even a single, long sightline shows a difference 
between the lower- z part and the higher-z part more than 
3cr (see Sect. 3.3 for more discussions). 

Similarly, £» c .i (13.6) derived from a larger number of 
lines is smaller, possibly approaching asymptotic values 
(Sample A: filled circles, filled squares and filled dia- 
monds). Both the & c i(13.6) values from Sample A and 
Sample B averaged (filled triangles) increase as z de- 
creases, although the uncertainties for Sample B averaged 
are rather large. Note again that the larger fluctuation 
of b Cj i(13.6) at < z > = 3.3 and 3.8 (keep in mind the 
smaller sample sizes at z > 3 at the same time). The 
smoothed b power-law fit also produces the large fluctua- 
tions of (r — l) s and 6 C , s (13.6) for the individual members 
of Sample B (see Table ||) . 

Our (r — I); and 6 Cj i(13.6) for Sample A are not con- 
sistent with the Schaye et al. values by more than 3cr. Our 
(r — 1); and 6 Cj i (13.6) for Sample B, however, are com- 
pletely in agreement with theirs. This result shows that 
deriving (T— l)i and 6 Cj i (13.6) depends largely on the num- 
ber of lines in the fit and is subject to the cosmic variance. 
Experiments with the fitted line parameters show that in 
general more than 200 lines in the fit stabilizes the results 
(cf. Schaye et al. 1999). The large fluctuation found from 
sightline to sightline in similar redshift ranges using the 
number of lines smaller than 200 could explain in part 
the previous contradictory conclusions on the evolution of 



(r — 1) and 6 C (13.6), which were usually based on a small 
number of lines and of sightlines. 

There are several known voids (or regions devoid of 
absorption lines with JVhi > 10 13 5 cm~ 2 ) in the present 
data. One void (z ~ 3.2) towards Q0302-003, two (z ~ 3.1 
and 3.3) towards Q0055-269 and three (z ~ 1.9, 2.2 and 
2.3) towards HE2217-2818 have been identified. The re- 
gions containing voids of Q0302-003 and Q0055-269 show 
a flatter (r — 1); by more than la (two lower squares 
at z ~ 3.1 in the upper panel of Fig. ||) and a higher 
6 Ci i(13.6) by more than 3cr (two higher squares at z ~ 3.1 
in the lower panel of Fig. ^ compared to the ones de- 
rived from Q0055-269 without the voids. The forest of 
HE2217-2818, however, does not show any significant dif- 
ference between the regions with the voids and without 
the voids. Similarly, the HE1 122-1648 forest without any 
recognizable voids does not show any significant differ- 
ence between the lower- z part and the higher-z part of 
the spectrum. 

Even though voids are produced by the enhanced radi- 
ations from local sources, the increase of the Hi photoion- 
izing photons in voids does not increase significantly the 
temp erature of the forest as whole (Haehnelt & Steinmetz 
1998 ). A fixed A^hi, however, corresponds to a higher over- 
density due to the increased photoionization. From the 
equation of state, a higher overdensity means a higher 
temperature. Therefore, if voids are produced by enhanced 
ionizations, there should be an increase of fo c .i(13.6) and 
a decrease of (T — 1); compared to the forest without any 
extra heating sources at similar z (Haehnelt & Steinmetz 
1998| Schaye et al. |1999| ). 

The lack of a significant difference in (T — 1); and 
& c a(13.6) from the HE2217-2818 forest at z ~ 2.1 sug- 
gests that the origin of the voids towards HE22 17-28 18 is 
differ ent. T his could be due to a density fluctuation (Heap 
et al. 200C) and/or due to a shock heating by gal actic i n- 
falls or by galactic winds (Theuns, Mo & Schaye 2001a| ). 

Simulations often show that some Lya lines are broad- 



ened by shock heating (Bryan & Machacek 200C; Theuns 
et al. 2001a). They also show that shock heating is not 



very important at z > 2 for the lower column density for- 
est considered here. In addition, the cutoff b method is 
less sensitive than other methods to shock-heated gas (cf. 
Theuns et al. |2001bp. 



One of the other main candidates for the fluctuations 
of (r — 1); and b c ,j (13.6 ) is t he Hen re ioniza tion at z ~ 3 
(Hae hnelt & Steinmetz 199S ; Songaila 1998 ; Schaye et al. 
1999|) . 



3.3. The Hen reionization 

There have been various claims on a possible Hen reion- 
ization at z ~ 3 from obs ervat ions (Rcime rs et al. 1997 



Songaila |1998[ Heap et al. |200C| ; Kriss et al. |200l|) . One of 

the consequences of the He II reionization is an increase of 
the t emperature of the Lya forest (Haehnelt & Steinmetz 
1998 ). In fact, various theoretical models have predicted a 



6 



Tae-Sun Kim, Stefano Cristiani and Sandro D'Odorico: The evolution of the physical state of the IGM 



(7 - 1) - (thus a flatter (T - 1)) and a higher T (thus 
a higher 6 C (13.6)) when the He ill bubbles surrounding 
the ionizing sources overlap and increase the radiations in 
the IGM (Ha ehnel t & Steinmetz |l99S| ; Ricotti et al. |200C ; 
Schaye et al. |2000|). 



To study the effect on the Lya forest caused by the 
He 11 reionization, we re-grouped the individual mem- 
bers of Sample B. For each z bin, individual members 
were divided into two subgroups according to their red- 
shift ranges, i.e. a lower-z part and a higher-z part. This 
method of grouping has been taken since the absorption 
lines from the same sightline might be correlated. The 
number of lines at each subgroup is sampled to be similar 
(see Table §. 

There is no significant difference in (r — l)i and 
6 Cj ;(13.6) between z ~ 2.0 and z ~ 2.1. There is, however, 
a significant difference by more than ~ 3er between the 
lower-z and the higher-z parts at z ~ 3.3 and at z ~ 3.8. 
The slope (r — l)i decreases at lower z for each bin at 
z > 3, while the 6 Ci ;(13.6) increases. In short, along each 
line of sight at z > 3, (r — l)i decreases and 6 Cj i (13.6) 
increases at z decreases, as expected from the He 11 reion- 
ization at z ~ 3. This trend, however, does not continue 
at z < 3, i.e. the HE1 122-1648 forest and the HE2217- 
2818 forest do not show the similar behavior. Thcuns ct 



al. (2001b) find a similar behavior of b values from the 
wavelet analysis using the same data presented here. This 
is interpreted as a result of Hell reionization at z ~ 3.3. 
They do find, however, a cold region, i.e. a region with 
a lower average b values than the adjacent regions along 
HE2217-2818. 

Deriving the cutoff Doppler parameters using a small 
number of lines introduces a large scatter (Fig. |J). 
Therefore, a flatter (r — l)j and a higher 6 c ,i (13.6) shown 
at the lower-z part of the bins at z > 3 (z ~ 3.1 and 3.6) 
might be an artifact of the small number of lines in the fit. 
More available lines in the fit, however, tends to increase 
(r — 1); and to decrease 6 C) i(13.6) as seen in Sect. 3.2, 
contrary to this result. 

The -/Vhi-6 distributions at z ~ 3.1 and z ~ 3.6 show 
a lack of lines with Ahi < 10 13 7 cm -2 and b < 20 km 
s _1 with respect to the ones at z ~ 3.4 and 3.9, while 
the Nm-b distributions for higher- A^hi lines are similar 
(the figures not shown). This lack of lower- N-ri and lower- 
b lines results in a flatter (r — 1); and a higher 6 c i (13.6) 
at z ~ 3.1 and 3.6. This, however, can not be caused by 
incompleteness due to line blending. The number of lower- 
JVhi lines are similar at z ~ 3.1 and z ~ 3.4, and at z ~ 3.6 
and z ~ 3.9, i.e. b values corresponding to lower- AThi lines 
are higher at z ~ 3.1 and z ~ 3.6 than those at z ~ 3.4 
and z ~ 3.9. In addition, line blending is more severe at 
higher z and should have resulted in a lack of lower- TVhi 
and lower-6 lines at z ~ 3.4 and z ~ 3.9 rather than at 
z ~ 3.1 and z ~ 3.6 as observed. What is observed at 
z ~ 3.1 and z ~ 3.6 shows the opposite behavior expected 
from line blending. In fact, the lower- A^hi and higher-6 
lines seen at z ~ 3.1 and at z ~ 3.6 are expected from 
the Hen reionization since it has a greater effect on the 



lower- ./Vh 1 forest than on the higher- A^hi forest (Theuns 
2001, private communication). 

Songaila (1998) found that there is an abrupt, sharp 
discontinuity in the ratio of Si iv column density to C IV 
column density, Agi IV /Aciv, at z ~ 3, below which 
Nsiw/Ncw is always less than ~ 0.07. The observed 
-Asiiv/Aciv by Songaila (1998) suggests that the UV back- 



ground is softer than a QSO-dominated background at 
z > 3 and becomes harder as expected from a QSO- 
dominated background at z < 3. This observation has 
been interpreted as the complete He 11 reionization by 
z ~ 3, i.e. the overlap of the He in bubbles. 

Fig. [| shows ATsirv/A^Civ of the Lya forest at iVm = 



10 



14-17 



cm 2 as a function of z from the QSOs in Table @ 
except from Q0000-263 (Kim et al. 2001, in preparation). 
We only include the H 1 systems having higher Lyman 
lines other than Lya, such as Ly/3, Ly7, etc. This selec- 
tion enables us to estimate an accurate H 1 column density 
and to assign Siiv and Civ to a Hi line more reliably. 
Note that there are no data points at 2.6 < z < 2.9 
at which the discontinuity of Ngi lv / Nqiv has been re- 
ported to be the largest. Although the bulk of the forest 
shows lower A^iiv/A/civ at z < 2.6, there are the forest 
lines with Asiiv/Aciv larger than 0.07. The similar re- 
sults have been reported at z > 2 by Boksenberg, Sargent 
& Rauch 



Bear in mind that no specific H I col- 
umn density range used in their figure 1 is given in their 
paper. In addition, note that their results are based on 
the individual, fitted components, while ours are based 
on the integrated absorption lines. Using the component- 
by-component analysis introduces an additional scatter 
in the diagram due to the different velocity structures in 
Si iv and Civ. The number of components for one inter- 
grated system in our data is in general from 1 to 3. The 
component-by-component analysis on our data also shows 
the similar trend found in Fig. |] from the integrated pro- 
files. 

Our observations do not suggest any abrupt change in 
Asiiv/A'civ at 1.6 < z < 3.6, i.e. no abrupt change in the 
softness parameter in general. This result suggests that 
some forest clouds are exposed to a soft UV background 
at z < 2.6. The lack of a sharp change in Nsiw/Ncw at 
z ~ 3 does not mean, however, that the He 11 reioniza- 
tion did not occur at z ~ 3. Rather, our Asiiv/Aciv indi- 
cates that ATgirv/A^civ might not be a good observational 
tool to probe the Hen reionization. What it probes is a 
softness parameter of the UV background. Increasing evi- 
dences of the contribution to the UV background from lo- 
cal, high-z galaxies would explain our high Nsuv/Ncrv at 
z < 2. 6 wit hout any difficulty (Giroux & Shull |l997j ; Shull 
et al. 119991; 1999 Bianchi, Cristiani & Kim 2001; Steidel, 



Pettini & Adelberger 2001). In addition, the lack of strong 



z-evolution of (r — 1) found for Ahi 



10 



13.7-14.5 



at z > 3 is in agreement with the lack of strong evolu- 
tion of N siIV /N C w at Ah 1 = 10 14 ~ 17 cm' 2 In s hort, the 
He n rei onization at z ~ 3 (Reimers et al. 1997 Heap et 
al. 2000 ; Smette et al. 2000 ) shows its impact on the Lya 



forest mainly at A^hi < 10 cm and might be very 
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Fig. 4. The ratio of Nsiw/Ncw as a function of z 
Only the Lya forest with A^hi < 10 
eluded. Filled circles represent a forest with Nq 

10 12.7-14 cm -2 

N Clv > 10 14 cm 
or lower limits. 



(in while open circles represent a forest with 
Open squares represent upper limits 



inhomogeneous. Its strength might not be as strong as 
previously suggested and Nsi lv / Nqw might not be a best 
observational tool to probe the He n reionization. 



forest from the He II reionization than their simulations 
suggest. When To is converted from the same simulations 
(their figure 2), To decrease as z decreases. 

For (7 — 1), we assume the conversion law between ]Vhi 
and 5 by Schaye (2001) assuming Jhm and T ~ 59.2 6 2 for 
thermally broadened lines. The conversion law is defined 

by 



Hi 



2.7 x 10 13 (1 + 5) 



1.5-0.26 (7-1) rp-0.26 
± r ■ 



l + Z 



4.5 



n h h 2 
0.02 



1.5 



0.4 



0.5 



T 



Hl.12 



0.16 



(2) 



where To = Tq^ x 10 4 K, the Hi photoionization rate 
Thi = Thi,i2 x 10 12 s -1 , fib is the baryon density, h is the 
Hubble constant divided by 10 0, and / g is the fraction of 
the mass in gas ( Schay e |2001 ) . We read To from figure 3 
of Schaye et al. ( 200C| ) and Thi from figure 8 of Haardt 
& Madau ( 1996Q , while we assume fib h 2 = 0.02 and / g = 
0.16. At < z> = 2.1, 3.3 and 3.8, (7 - 1) - 0.417, 0.386 
and 0.441. Assuming Jhmj there is no clear z-evolution of 
(7 — 1) within large uncertainties. Our values and Schaye 
et al.'s ((7 - 1) - 0.4, 0.35 and 0.25 at < z > = 2.1, 3.3 
and 3.8 assuming Jhm) agree at < z > = 2.1, 3.3, while 
our value is a factor of 1.8 larger than theirs at < z >= 
3.8. Our (7 — 1) values agree with those of McDonald et 
al. at the similar z ranges within uncertainties, although 
their simulations do not assume Jhm. 



4. Discussion 



5. Conclusions 



Converting (r — 1) and 6 C (13.6) to the corresponding 
(7 — 1) and To depends on many uncertain parameters, 
such as the UV background and the reionization history. 
If the UV background is dominated by QSOs without 
any extra heating, Tq decreases and (7 — 1) increases as 
z decreases, until they approach asymptotic values. The 
6 C (13.6) value, however, still increases with decreasing z, 
until it approaches an asymptotic value and finally de- 
creases again at z < 2 (Schaye et al. 1999). 

Instead of converting to To, we compare o ur 6 c , j(13.6) 
with figure 3 and figure 4 of Schaye et al. ( |200C ) from 
their simulations assuming the QSO-dominated Haardt- 
Madau UV background, J H m (Haardt & Madau |1996| )p|. 
Our 6 C (13.6) is in agreement with their simulated 6 C (13.6) 
values without extra He 11 heating, although the difference 
increases as z decreases (ours being a factor of 1.1 lower 
at z ~ 2). It is not clear what causes their numerical 
simulations without the Hell reionization produce their 
6 C (13.6) similar to ours which suggest the extra heating. 
One of the explanations could be a weaker effect on the 



The previous section suggests there is a contribution from 
sources other than QSOs to the UV background. Estimates 
of contributions from galaxies depend on many factors, such 
as the galaxy luminosity function and the galaxy e missiv ities. 
Including a contribution from galaxies, Shull et al. ( 199S ) de- 
rive a UV background decreasing from 2 ~ 2.5 to 2 ~ 1, similar 
t o the QSO-dominated Jhm - On the other hand, Bianchi et al. 
(2001) find a non-decreasing J from z ~ 5 to z ~ 1. 



Using the new, large dataset from high S/N, high reso- 
lution VLT/UVES data combined with one Keck/HIRES 
QSO in the literature, the minimum cutoff Doppler pa- 
rameter as a function of Nni, b c (Nni), of the Lya forest 
has been derived at < z > — 2.1, 3.3 and 3.8. We have 
found: 

1. When incompleteness is accounted for, the derived 
(r — 1) is consistent with no-z evolution with a sugges- 
tion of lower (r — 1) at z ~ 3.1, while 6 C (13.6) increases 
as z decreases. These results suggest that (7 — 1) shows 
no z-evolution within uncertainties and that To de- 
creases as z decreases, assuming the QSO-dominated 
UV background from Haardt & Madau ( |1996[ ). 

2. The (r — 1) and 6 C (13.6) values show a large fluctua- 
tion when derived using a subsample even at the simi- 
lar redshifts. The fluctuation is larger at z > 3 than at 
z < 3. Although smaller number of lines in the fit could 
result in this fluctuation, we use a similar number of 
lines to derive (r — 1) and 6 C (13.6) for each subsample. 
This result suggests that there might be a large fluctu- 
ation in the IGM temperature along different sightlines 
even at similar z. 

3. There is a suggestion of a flatter (r — 1) and higher 
fe c (13.6) along each line of sight as z decreases at z > 3 
more than 3cr, probably due to the He 11 reionization. 
At z ~ 2.1, there is no such a significant trend along 
each line of sight. This result occurs due to the lack of 
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lower- A^hi and smaller-^ lines in the lower- z part of the 
spectra at z > 3. Our result implies that the impact 
of the Hen rcionization on the T.yrv forest, might he 



Shull, J. M., Roberts, D., Giroux, M., Penton, S. V., Fardal, 
M. A., 1999, AJ, 118, 1450 



Smette, A., He ap, S., Williger, G. M. et al, 2000, [astro- 



malnly on the lower- /V Hl forest and that its significance 



ph/0012193| 



might be smaller than previously suggested. 

The lack of strong discontinuity of Nsiw / New at 

Nm = 10 14 " 17 cm- 2 at z ~ 3 suggests that N SUv /N C w 



might not be a good observational tool to probe the 
He II reionization and that the UV background might 
be strongly contributed by local, high-z galaxies at 
z < 3. 
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